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Abstract. This article discusses recent attempts to provide a deeper understanding of the thermoreversible 
" gel" state of colloidal matter and to unravel the analogies between gels at the colloidal level and gels at the 
molecular level, commonly known as network-forming strong liquids. The connection between gel-forming 
patchy colloids and strong liquids is provided by the limited valence of the inter-particle interactions, i.e. 
by the presence of a limit in the number of bonded nearest neighbors. 

PACS. 82.70.Gg Gels and sols - 82.70.Dd Colloids 61.20.Lc Structure of liquids; time-dependent proper- 
ties; relaxation 

1 Introduction the experimental observation time. Despite the apparent 

simple definition, the essence of the (colloidal) gel state 
Colloidal gels are an arrested state of matter at low den- ig gtm under ag connrmed by the pub i ica ti n of 
sity or packing fraction, i.e. an arrested state in which geyeral recent reyiewf , summarizing the work which has 
a negligible fraction of the total volume is occupied by been donc in thfi lagt vears |l | 2 | 3 l 4 l 5] . 
particles. Shear waves of large wavelength can be trans- 
mitted trough the gel via the bonded-particle percolat- In this article j wiU mostly focus on the expected be- 
ing structure. In the gel, particles are tightly bonded to havior of recently newly synt hesized colloids: nano and 
each other. This requires a thermal energy significantly micr0SC0 pi c particles with a chemically modulated sur- 
smallcr than the bonding energy, so that the lifetime of fac6; to provide colloids with valence [SEH] . This rel- 
the thermo-reversible bonds is comparable or longer than evant synthesis effort aims to gene rate super atoms - 
Send offprint requests to: atoms at the nano and micro-scopic level — to reproduce 



2 Francesco Sciortino: Gel-forming patchy colloids and network glass formers: Thermodynamic and Dynamic analogies 

and extent the atomic and molecular behavior on larger A series of recent studies have focused on the role of the 
length scale, moving from the colloid/hard- sphere-model valence in determining the field of (meta)stability of the 
analogy to more complex cases. Here I exploit one of these fluid disordered phase [TT|12I13| . The field of (meta)stability 
analogies: (thermoreversiblc) gel-forming patchy colloidal can be defined as the region in the temperature-density 
particles and associated liquid, suggesting that network plane in which the equilibrium homogeneous fluid (or liq- 
strong glass forming liquids are the atomic and molecular uid) phase can exists for a time sufficiently long to allow 
counterpart of patchy thermoreversible colloidal gels. for an experimental (or numerical) observation. Even if 

crystallization can be neglected (due to the presence of 
2 The Role Of the Valence a nucleation barrier which can sometime be longer than 

the experimental observation time), the field of stability of 



The connection between patchy colloidal gels and network 
glass-forming liquids is provided by the limited valence 
of the inter-particle interactions, i.e. by the presence of a 
limit in the number of bonded nearest neighbors. The lim- 
ited valence can be brought in or via strong directional in- 
teractions (hydrogen bond is the classical molecular exam- 
ple), via interplay of isotropic competing interactions [10 
or due to specific three-body interactions (as in silicon). 
At odd with standard atomic or molecular liquids, in lim- 
ited valence systems packing is not the driving force deter- 
mining the structure of the system. The number of near- 
est neighbors in the liquid state is significantly smaller 
than twelve and the hard-sphere model does not play the 
role of reference system. Crystal structures are character- 
ized by solids with low coordination number, as in the 
four-coordinated hexagonal ice (H2O) and quartz (S1O2) 
cases. These crystal structures can be considered empty. 
The packing fraction of touching hard-spheres arranged 
in the diamond structure is indeed only w 0.34, less than 
half of the packing of the hard-sphere face-centered cubic 
structure. 



the fluid phase is still limited by the glass-line (preventing 
equilibration of the sample on the experimental time scale) 
and by the liquid-gas effective spinodal (preventing the 
equilibration of the sample in an homogeneous structure). 
In the case of particles interacting with spherical symmet- 
ric potentials (e.g. simple liquids) or non-symmetric but 
not valence-limited potentials, the glass line intersects the 
gas-liquid spinodal on its right side (the liquid side) [2], 
as schematically shown in Fig. [TJ-(top). Thus, at low tem- 
perature it is not possible to generated disordered homo- 
geneous (as opposed to phase separated) arrested states. 
Phase separation covers a large region of the phase dia- 
gram. For square-well type interactions (for which a clear- 
cut estimate of the packing fraction <j) can be provided 
since the hard-core distance is properly defined) the crit- 
ical packing <f> c moves from <j) — 0.13 in the infinite in- 
teraction range limit (Kac model) up to <f> — 0.27 in the 
infinitesimal interaction range limit (Baxter model) [TS]. 
Since the glass line meets the spinodal line at density 
larger than </> c , it follows that it is not possible to ob- 
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serve homogeneous arrested states at densities lower than in Fig. [T|-(bottom) [ll|12| . This has been demonstrated 
<f> c . theoretically (via application of the Wertheim thermody- 
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namic perturbation theory [16117] ) and via simulations for 
a series of simple patchy models, in which the interaction 
between particles is composed by an hard-core comple- 
mented by several patch-patch interactions, under the con- 
dition of single-bond per patch [TTl 12] . It has been shown 
that the width of the liquid-gas coexistence region pro- 
gressively reduces on decreasing the valence, disappear- 
ing when the average functionality approaches two. For 
small valence, the Wertheim theory has been shown to be 
very accurate, providing a parameter free description of 
the system in a wide density and temperature region. The 
significant shrinking of the gas-liquid region is particularly 
important in the present context. Indeed, it becomes now 
possible to reach low temperatures without encountering 
phase separation in conditions such that packing is not an 
effective arrest mechanism. This region of low tempera- 



Fie. 1. Schematic phase diagram of the disordered phases for , /. £ , , tjuj- \ j-ci j 
& 1 6 v ture (i.e. of strong long-lived bonding) and of low densi- 

the case of particles interacting with non- limited (top) and lim- ,. . . .. . jN 

ties (i.e. m condition where packing is not relevant) can 

not be accessed in spherical isotropic potentials since it is 
pre-empted by the phase-separation process. In the case 
of limited valence particles instead, even at intermediate 
densities, at low T most of the particles are fully bonded 
and there is no energetic driving force toward a further 
increase in the local density. The shrinking of the unsta- 
ble region depends on the average valence. In the limit of 
average valence approaching two (which can be realized 
for example by a binary mixture of particles with two and 
three patches, in the limit of vanishing concentration of 



ited (down) valence potentials. In the standard (top) case, the 
glass line hits the gas-liquid spinodal at large densities. In the 
limited valence case, the shrinking of the gas-liquid unstable 
region opens up a new region in which a stable extensively 
bonded network develops. In this region, the dynamics is the 
one of strong glass-formers. 

Reduction of the valence provides an effective mecha- 
nism to shrink the gas- liquid coexistence region, opening a 
large region of intermediate densities in which there is no 
driving force for phase separation, as schematically shown 
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three-functionalized particles) the region of instability ap- very often very well defined interactions. With minor am- 
proaches zero. For average valence approaching two, a sta- biguities, it is possible to define the presence of a bond 
blc homogeneous empty liquid phase of vanishing density between two arbitrary particles (in simulations) or detect 
can be realized. The structure of the system is composed the fraction of bonded pairs with spectroscopic measure- 
by long chains of bi- functional particles cross-linked by the ments. The study of the connectivity properties of limited 
three- functional ones. The possibility of satisfying all pos- valence systems shows the presence of a bond percolation 
sible bonds in a wide region of densities provides thermo- line, passing always slightly above the critical point. At 
dynamic stability to the structure. A detailed discussion small valence, the Cy-max line lies above the percolation 
of this limiting case can be found in Ref. [T5] . line [23 18, The opposite behavior is observed at large va- 
lence 

Recently it has been reported that universality in the 

behavior of different limited- valence models can be re- T , j- r • i j i -j.i 

1 now turn to discuss properties of simple models with 

covered using the valence to discriminate between differ- , r , , , .,, 

° valence four, to make contact with the most common net- 

ent behaviors. More precisely, when the leading interac- ,„. . . „ ,., , . ., 

° work forming glass formers, beveral simple patchy models 

tion is provided by bonding (i.e. when the residual in- r r . , . . , -, , . ,. , 

^ j o \ j j- our coordinated particles have been recently studied, 

teraction is essentially condensed in a hard-core repul- .,, . r n ,. , . , e , , , , 

J with the aim of finding out the universal features related 

sion) the location of the critical point is controlled by . , mi ii-ijj.it-j.ji 

' r J to the valence, ihese models include the limited valence 

a (valence-dependence) critical value of the second virial ln j wat mttftti ji • •,• i i r 

v ^ ; square-well model (N max ) [11 24, the primitive model for 

coefficient^. This scaling generalizes the Noro-Frenkel water (PM W) ESMOZl , the primitive model for silica (PMS) 

law of corresponding states for spherical interactions HH]. the floating bond modeim the modd for four _ 

The scaling is associated to the possibility of expand- , mT , , , . i aa i aj . . , , ,. c 

° t- t- armed DJNA dendnmcrs[30 31J. All these models satisfy 

ing the partition function in term of pair-wise bond free - , , , , , ^ c ■ j-v. j- j. 

° the single- bond per patch condition. Defining the distance 

energy [22]. between a pair of bonded particles as unit of length and 

The thermodynamic of systems with limited valence is the critical temperature as unit of temperature, all these 

characterized by the presence of a temperature at which models are characterized by a similar phase diagram (Fig. [21. 

the constant volume specific heat Cy has a maximum. The The gas-liquid unstable region is limited to scaled densi- 

locus of Cy extrema in the phase diagram is a clear ex- ties lower than 0.6. Above that density, the liquid is 

perimentally detectable evidence of the leading role of the stable and the structure of the liquid is characterized by 

bonding interaction in the thermodynamic of the system, an extended network of four coordinated particles (being 

since it can be interpreted as arising from the equilibrium the liquid state well within the percolation region). In all 

between broken and formed bonds. Moreover, bonds are the single-bond per patch models, a well defined ground 
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state exists, provided by the geometrically fully bonded 
configuration. In all the studied models, the fully bonded 
state can be accessed in a region of densities whose width 
is controlled by the angular constraints [32 . The approach 
to the ground state energy along isochoric paths is well de- 
scribed by an Arrhenius law in temperature. Universal fea- 
tures also control the particle dynamics in these systems. 
At low T, the T-dependence of the diffusion coefficient D 
and of the viscosity r\ approaches an Arrhenius law, the 
characteristic behavior of so-called strong glass-formers. 

An important consequence of the possibility of ap- 
proaching continuously the ground state of the system 
(which in these models is a-priory known, being controlled 
only by the number of patches per particles and by the 
single-bond per patch condition) provides a proof that 
there is no vanishing of the configurational entropy at any 
finite temperature |33I3^] . at least in this class of models. 
More specifically, there is not finite Kauzmann temper- 
ature. Even more interestingly, in the inherent-structure 
thermodynamic formalism [35 36 , the configurational en- 
tropy (which counts the number of distinct configurations 
with the same energy) can be evaluated numerically with- 
out any approximation [33 . Indeed the number of distinct 
configurations can be formally associated to the number of 
topological distinct bonding patterns with the same num- 
ber of bonds and the vibrational entropy to the phase- 
space explored within each bonding pattern. Results sug- 
gest that the ground state is degenerate and that the de- 
generacy is a function of the geometric constraints limiting 
the bonding angles. The degeneracy of the ground state 
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Fig. 2. Gas- liquid spinodal line and (lowest available) iso- 
diffusivity line for five different models of patchy particles, in 
the scaled variables T/T c and scaled p a caied- In all models, 
the gas-liquid unstable region extends for p sca ied < 0.6. The 
iso-diffusivity lines, providing an estimate of the shape of the 
glass line, are all mostly parallel to the T-axis, bending up only 
on approaching large densities. Scaled densities are defined as 
number of particles divided by the volume, measured in units of 
the average distance between two bonded particles. Data are 
from: (PMW) [26], (PMS) [32], (AUxx)[EEl24], floating bond 
model [29 , the model for four-armed DNA dendrimers 31:. 

reflect the number of distinct fully bonded configurations, 
differing in the bonded-ring statistics. 

The class of valence-four systems includes two of the 
most abundant liquids on earth, water and silica. For these 
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materials, several continuous potentials have been devel- 
oped and extensively studied, providing a reasonable re- 
production of several thermodynamic and dynamic prop- 
erties of the real substances. These classical potentials, 
based on electrostatic and Lennard-Jones type of inter- 
actions, are able to reproduce the tetrahedral structure 
characteristic of water and silica under the appropriate 
density and temperature conditions. At low T, the diffu- 
sion coefficient shows an Arrhenius T dependence. It is 
thus worth to compare the behavior of these more compli- 
cated and accurate potentials with the generic features 
observed previously in the simpler one-bond per patch 
models. For all models, the generic features observed in 
the patchy colloidal particles models are confirmed. Even 
the structure of the tetrahedral network generated by the 
ionic model for silica is astonishingly similar to the one 
generated by the short-ranged patchy model for silica[32 . 
The gas-liquid unstable region is located in the same re- 
gion of scaled densities and the iso-diffusivity lines are 
essentially controlled by temperature (parallel to the T- 
axis), as shown in Fig. [3] 

3 Conclusions 

The reviewed works point to the relevant role of the lim- 
ited valence in the formation of a gel state of matter. 
Dynamic arrest at low densities in an homogeneous sys- 
tem requires a very strong inter-particle attraction (signif- 
icantly larger than the thermal energy ksT). Under these 
conditions, phase separation can be beaten only when 
minimization of the energy does not require an increase in 
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Fig. 3. Phase diagram, including iso-diffusivity lines, of three 
different models: (a) the BKS [37] model for silica, (b) the 
ST2[38] and (c) the SPC/E[39] models for water. Dashed lines 
are the gas-liquid spinodal, full lines are the slowest available 
iso-diffusivity loci. Different isodiffusivities lines are spaced by 
one order of magnitude in the value of the diffusion coefficient. 
The diagrams are rebuilt using published data. For the BKS 
potentials data are from Ref. [40141] . For the ST2 potentials 
data are from Ref. [42] and for the SPC/E potential data are 
from Ref. 43_. In the case of water, distances are measured in 
units of 2.8 A (the HB distance), so that 1.0 g/cm 3 correspond 
to a scaled density of 0.734. In the case of silica, distances are 
measured in units of (3.1 A) (a characteristic value for the Si-Si 
distance), so that 1.0 g/cm 3 correspond to a scaled density of 
0.3. 
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the local density. Interacting patches which can not form models the ground state is continuously approached on 

multiple bonds are particularly effective in this respect, cooling, proving the absence of any thermodynamic tran- 

In this spirit, it is possible to rationalize the behavior sition associated to the dynamic arrest. No finite Kauz- 

of polymeric systems with a small fraction of functional- mann temperature exists. 



ized monomers, the behavior of particles interacting with 
lock-and-key type of interaction, very often encountered 
in biological self- assembly [44 45 46 ) as well as associated 
liquids and hydrogen bonded compounds [47 . 



In fragile glass- formers, homogeneous arrested states 
are observed only at very high densities, when excluded 
volume caging becomes dominant. At low density, phase 
separation prevents the possibility of homogeneous arrest. 
It is thus tempting to speculate that dynamics can be 
mostly interpreted in term of hard-sphere dynamics, simi- 
larly to the role played by the hard-sphere model as refer- 
ence system for structural and thermodynamic quantities. 
In this picture, energy and temperature play a role mainly 
via the change in the parameters of the reference state. If 
this is the case, fragile and strong liquids are two distinct 
(limiting) classes of materials. 



There are several questions arising from the analogy 
developed between network forming liquids and patchy 
particles. In particular I like to call attention on the fact 
that the reduction of the valence appears to be critical 
in opening up an intermediate region of densities where 
bonding (as opposed to packing) becomes the leading in- 
teraction. The slowing down of the dynamics in this re- 
gion is found to be Arrhenius and hence typical of the 
so-called strong glass formers. Arrest is driven by the for- 
mation of a network of long-lived bonds, i.e. by the for- A final comment regards the behavior of liquid water, 
mation of energetic cages, as opposed to the more famil- As discussed, the gas-liquid unstable region is analogous to 
iar excluded volume cages characteristic of fragile glasses, the one observed in tetrahedral coordinated patchy mod- 
In the studied models, the slowing down of the dynam- els. Interestingly, several water models are characterized in 
ics is indeed triggered by the formation of bonds and the region of intermediate densities by an additional region 
the Arrhenius dynamics arises when most of the possi- of instability (the second liquid- liquid critical point 48J) 
ble bonds are formed. The simplicity of the Arrhenius separating two disordered liquid structures with different 
dynamics appears to be connected to the elementary lo- densities. The low-density liquid has indeed the proper- 
cal independent process of bond-breaking. In the simple ties of an empty liquid structure while the high density 
model investigated [26 32 30 31 , the diffusion process is one appear to be characterized by a much larger density, 
controlled by the very small concentration of free par- It is not a coincidence, I believe, that dynamics in the 
tides (fully unbonded) diffusing in an empty matrix of low density liquid is Arrhenius (strong liquid) while it is 
fully bonded (and hence not reactive) particles. In these super-Arrhenius (fragile) in the high density one [49 . 
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